Adsorption and dissociation of hydrogen molecules on bare and functionalized carbon 

nanotubes 
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We investigated interaction between hydrogen molecules and bare as well as functionalized single- 
wall carbon nanotubes (SWNT) using first-principles plane wave method. We found that the binding 
energy of the H2 physisorbed on the bare SWNT is very weak, and can be enhanced neither by 
increasing the curvature of the surface through radial deformation, nor by the coadsorption of Li 
atom that makes the semiconducting tube metallic. Though the bonding is strengthened upon 
adsorption directly to Li atom, yet its nature continues to be physisorption. However, the character 
of the bonding changes dramatically when SWNT is functionalized by the adsorption of Pt atom. 
Single H2 is chemisorbed to Pt atom on the SWNT either dissociatively or molecularly. If Pt-SWNT 
bond is weakened either by displacing Pt from bridge site to a specific position or by increasing 
number of the adsorbed H2, the dissociative adsorption of H2 is favored. For example, out of two 
adsorbed H2, first one can be adsorbed dissociatively, second one is chemisorbed molecularly. The 
nature of bonding is weak physisorption for the third adsorbed H2. Palladium also promotes the 
chemisorption of H2 with relatively smaller binding energy. Present results reveal the important 
effect of transition metal atom adsorbed on SWNT and advance our understanding of the molecular 
and dissociative adsorption of hydrogen for efficient hydrogen storage. 
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I. INTRODUCTION 

Fuel cells have been a real challenge for clean and ef- 
ficient source of energy in diverse fields of applications 
with different size and capacity range. Once hydrogen 
molecule is chosen as potential fuel, its storage, easy dis- 
charge for consumption and dissociation into hydrogen 
atoms to complete the current cycle in the fuel cell to 
produce the desired electromotive force involve several 
problems to be solved. Dillon et alX have pioneered the 
idea that carbon nanotubes can be efficient, cheap and 
rechargeable storage medium for small-scale fuel cells by 
estimating 5-10 weight percent (wp) H 2 adsorption in 
single- wall carbon nanotubes (SWNT). Later, Ye et alwk 
and Liu et aim- obtained H 2 storage capacities of 8.2 and 
4.2 wp, respectively. Unfortunately, recent studies fur- 
ther exploring this idea have come up with controversial 
conclusions^*SiLSiS. In the mean time, adsorption of al- 
kali atoms on SWNTs have been proposed to enhance 
the RVuptak o 9 ! 10 . Nevertheless, carbon nanotubes have 
high surface-volume ratio and their functionalization to 
render them feasible for hydrogen storage through cover- 
age of suitable adatoms has remained to be explored. 

In this papesii, we addressed following questions in 
order to clarify controversial issues related to the stor- 
age of hydrogen molecule on carbon nanotubes: Can 
H 2 molecule be adsorbed on the SWNT? What is the 
nature and strength of the bonding? Can the strength 
of the bonding be modified either by changing the cur- 
vature of the surface or by the coadsorption of metal 
atoms? In particular, can the H 2 -uptake on SWNTs be 
promoted by their functionalization through transition 
elements? How can H 2 molecule be dissociated? To an- 



swer all these questions we investigated the interaction 
between H 2 molecule and bare, radially deformed and 
foreign atom adsorbed SWNTs by carrying out calcu- 
lations within Density Functional Theory (DFT)A£. Our 
results not only advance our understanding of H 2 adsorp- 
tion on carbon nanotubes, but also suggest new ways for 
efficient hydrogen storage for rapidly growing research on 
fuel cells. 

Our calculations have been carried out using 
first-principles plane wave method and ultrasoft 
pseudopotentialsiii4 within Generalized Gradient Ap- 
proximation (GGA)a*. Adsorption and dissociation of 
H 2 is treated within the supercell geometry by opti- 
mizing atomic positions. The weak attractive Van der 
Waals (VdW) interaction becomes crucial in calculating 
binding energies of weak physisorption bonds, but is 
not well represented in DFT using GGA 16 . Therefore, 
in the case of physisorption, weak and attractive VdW 
interaction energy, EvdWi is obtained from the Slater- 
Kirkwood approximation™! using the asymptotic form 
of the Lifshitz's equation 18 . However, this approach can 
not be suitable to determine the contribution of VdW 
interaction in the chemisorption of molecules. In the 
present calculations we take the zigzag (8,0) SWNT as 
a prototype tube. 



II. ADSORPTION OF H2 ON BARE AND 
RADIALLY DEFORMED SWNT 

To clarify whether H 2 can form stable bonding on the 
outer or inner surface of a SWNT, we calculated the 
chemical interaction energy between H 2 and the outer 
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FIG. 1: Variation of chemical interaction energy Ec between 
SWNT and H2 molecule as a function of distance d between 
them. Two cases, namely adsorption to bare and radially 
deformed SWNT are shown by dashed and continuous lines, 
respectively. In calculating both curves, atomic structures 
corresponding to d — > 00 have been used without relaxation. 
Dash-dotted line indicate zero of chemical interaction energy. 
Optimized distance for two cases are indicated by arrows. 



surface of the (8,0) SWNT at different sites {i.e. H-site, 
above the hexagon; Z- and A-site above the zigzag and 
above the axial C-C bonds; T-site, a bridge site between 
two adjacent zigzag C-C bonds) as a function of spacing 
d. The chemical interaction energy is obtained from the 
expression, E c {d) = E T [H 2 + SWNT,d]-E T [SWNT}- 
Et [H 2 ] , in terms of the total energies of bare nanotube 
(E T [SWNT]), free H 2 (E T [H 2 }), and H 2 attached to 
SWNT at a distance d (E T [H 2 + SWNT, d]). Here E c < 
corresponds to an attractive interaction. The stable 
binding occurs at the minimum of Ec{d) + Ey dw{d) , the 
negative of it is denoted as the binding energy Ei } . The 
binding is exothermic when Ef, > 0. In Fig. ^ we show 
the variation of Ec{d) calculated for unrelaxed atomic 
structures at the H-site. Once the atomic structure of 
both SWNT and H 2 molecules are relaxed the mimimum 
value of Ec(d) is found to be -27 meV at do = 3.1 A at 
the H-site. Minimum values of Ec{d) calculated for A-, 
Z-, and T-sites are also very small and comparable to 
that of H-site. The long range VdW interaction energy 
calculated for the H-site at d Q is Eydw ~-30 meV. Then 
the binding energy associated with H 2 molecule adsorbed 
at H-site is calculated to be E\, ^57 meV. This is a small 
binding energy and indicates physisorption. 

Previous studies-^ revealed that the physisorption of 
individual H 2 molecules with Ef, > cannot occur on the 
inner wall of SWNT. Hydrogen molecules prefer to stay 
either at the center of the tube or to form some cylindrical 
shells inside depending on radius of the tube. Owing 
to increased H-H interaction at high coverage, even the 
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FIG. 2: Atomic configuration, energy band structure and 
local density of states (LDOS) calculated for the coadsorption 
of H2 molecule and single Li atom. Two cases correspond to Li 
atom chemisorbed on the external and internal surface of the 
(8,0) zigzag SWNT. Zero of energy is set at the Fermi level, 
Ef. LDOS calculated at Li and H2 are shown by continuous 
and dotted lines, respectively. Metallized SWNT bands are 
indicated by arrows. 



atomic hydrogen cannot form stable structure when it 
is adsorbed to the inner wall of small radius tube. For 
(8,0) tube we found that H 2 is trapped and stabilized at 
the center of the tube with a repulsive interaction energy 
Ec = +0.34 eV. The implementation of H 2 inside the 
tube having radius in the range of 3 A is expected to be 
hindered by this repulsive interaction. 

Earlier it has been shown that binding energy of for- 
eign atoms adsorbed on SWNT increase with increasing 
curvature^. Tada et al£L have argued that the poten- 
tial barrier associated with the dissociative adsorption of 
H 2 on SWNT is lowered with increasing curvature of the 
tube. It has been proposed that the potential barrier for 
the dissociation of H 2 adsorbed in the interstitial region 
between tubes can be lowered by applying radial defor- 
mation to the rope or SWNT 21 . It is known that under 
radial deformation the circular cross section changes and 
consequently the curvature varies at different locations 
on the surface. Motivated with those effects of curvature, 
we examine whether the attractive interaction energy Ec 
can be enhanced by changing the curvature of the tube 
via radial deformation. Radial deformation is realized 
by pressing the tube between two ends of a given diame- 
ter. It, in turn, changes the circular cross section of the 
bare tube with radius R to an elliptical one with ma- 
jor and minor axes 2a and 2b, respectively. The atomic 
structure of (8,0) tube is optimized under radial strain 
e r = (R — b)/R ~0.3 by fixing row of carbon atoms at 
the end of minor axis. The deformation is reversible so 
that the tube goes back to its original, undeformed form 
upon the release of radial strain^. The deformation en- 
ergy (that is the difference between the total energies of 
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deformed and undeformed SWNTs) is calculated to be 
En =1.4 eV per unit cell. We examined whether the 
binding energy of H2 molecule changes under the radial 
deformation of SWNT. Figure 2] shows the variation of 
Ec(d) for H2 approaching toward the high curvature site 
of the tube (i.e. one end of the major axis) at the H-sitc. 
The minimum value of the attractive interaction, Ec is 
-30 meV and occurs at d =2.9 A. For e r =0.25, we also 
obtained very small enhancement of Ec(d). Hence, the 
effect of curvature on the binding energy of H 2 is negli- 
gible due to relatively large d . This result also suggests 
that the physisorption energy does not vary significantly 
depending on the radius of SWNT. Apparently, the bind- 
ing of H2 on the outer surface of SWNT is weak and cor- 
responding physisorption energy is small. The binding 
cannot be enhanced by increasing the curvature locally 
through radial deformation. Curvature effect or radial 
deformation may be significant at small d when H2 is 
forced towards SWNT surface. 



III. COADSORPTION OF HYDROGEN 
MOLECULE AND LITHIUM ATOM ON SWNT 




(8,0)+Pt (8,0)+2Pt (8,0)+3Pt 



E=2.92eV/Pt E =2.84 eV / Pt E,=3.42eV/Pt 

b l> I' 

FIG. 3: (a) Atomic configurations for single, double and triple 
Pt atoms adsorbed on the (8,0) SWNT. Average binding en- 
ergy of adsorbed Pt atoms Eb and bond distances are indi- 
cated. 

d decreased to 2.1 A. Briefly, the coadsorption of Li does 
not have any indirect effect on the binding of H2 , but the 
energy associated with direct binding to Li is enhanced. 
However, the nature of bonding remains physisorption in 
direct and indirect cases. 



Next, we examine whether the binding of H 2 is en- 
hanced by the coadsorbed foreign atoms. To this end we 
first consider Li atom adsorbed on the (8,0) SWNT, since 
the adsorption of an alkali atom has been proposed to en- 
hance the H2-uptake 9 ' 10 . Li atom is chemisorbed at the 
H-site, 1.5 A above the surface of SWNT with a binding 
energy of 0.8 eV. Self-consistently calculated electronic 
structure shown in Fig. [2] reveals that chemisorbed Li 
atoms donate their 2s-valence electrons to the lowest con- 
duction 7r*-band so that the semiconducting (8,0) SWNT 
(having band gap E g =0.6 eV) becomes metallic. This 
is a behavior common to the other alkali atoms adsorbed 
on SWNTs^. In order to examine the indirect effect of 
coadsorbed Li we considered H2 as attached at the oppo- 
site site to Li. The optimized structure of the physisorbed 
H2 is shown in Fig. [21 together with relevant structural 
parameters. We found Ec has a minimum value of -35 
meV at d a =3.4 A. Similar study has been also performed 
for Li atom adsorbed on the inner wall of SWNT while 
H2 is on the external wall directly above the coadsorbed 
Li as shown in Fig. [21 In this adsorption configuration 
minimum value of Ec practically did not change. The lo- 
cal density of states calculated on Li atom and H 2 refuse 
the possibility of any significant interaction between ad- 
sorbates. As a result, our calculations for both external 
and internal adsorption of Li rule out any indirect effect 
of coadsorbed Li to enhance the binding of H2 on SWNT. 
The occupation of empty conduction band by the alkali 
electrons and hence metallization of SWNT did not affect 
the bonding of H2. These results are in agreement with 
the first principles calculations by Lee et alw&. However, 
the effect of Li on the adsorption of H2, whereby H2 is 
attached directly to Li atom is found significant. The 
minimum value of Ec has increased to -175 meV, while 



IV. COADSORPTION OF HYDROGEN 
MOLECULE AND PLATINUM ATOM ON SWNT 

A single transition metal atom adsorbed on the outer 
surface of SWNT has shown interesting properties, such 
as high binding energy and magnetic ground state with 
high net magnetic moment. For example, transition ele- 
ment atoms (Ti, V, Cr, Mn, Fe, Co, Pt, etc.) have crucial 
adsorption states on nanotubes2£ and some of them (Ti, 
Ni, Pd) form continuous or quasi-continuous metal coat- 
ing on the SWNT212£. As for Pt atom, it is known to 
be a good cathalist in various chemical processes. While 
SWNTs offer high surface/volume ratio, the interaction 
between H2 and Pt atom adsorbed on SWNT may be of 
interest. Now we investigate the character of the bonding 
between H2 and Pt adsorbed on SWNT and address the 
question of how many H2 molecules can be attached to 
an adsorbed Pt atom and how strong is the binding. 



A. Adsorption of Pt atoms on SWNT 

We first examine the adsorption of Pt atom(s) on (8,0) 
SWNT. The character of the bonding has been investi- 
gated by placing Pt atoms on the A-sites of (8,0) tube 
(that is known to yield highest binding energy^) and 
then by optimizing the structure. Three different ad- 
sorption configuration have been examined, namely one, 
two and three Pt atoms adsorbed on the adjacent sites to 
represent a small cluster on SWNT as described in Fig. [3] 
Calculated binding energies of Pt atoms have increased 
as the number of Pt atoms increases from one to three 
in the same neighborhood. On the other hand, the C- 
Pt distance gradually increases with increasing number 
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FIG. 4: Optimized binding configuration of H2 molecules ad- 
sorbed to a free Pt atom, (a) Dissociative adsorption of a sin- 
gle H2 molecule, (b) The first H2 is dissociatively, second H2 
molecularly adsorbed, (c) Two H2 are molecularly adsorbed. 

(d) Two H2 are molecularly, one H2 dissociatively adsorbed 

(e) Four different configurations related with the adsorption 
of two H2 to the same free Pt atom and corresponding four 
minima on the Born-Oppenheimer surface described schemat- 
ically. Binding energy of each adsorbed additional hydrogen 
molecule Et,, total energy with respect to constituent atoms 
Et and bond distances are indicated. 



of Pt atoms adsorbed in the same neighborhood. This 
paradoxical situation can be understood by the increas- 
ing Pt-Pt coupling, which happens to be comparable to 
C-Pt coupling derived by the Pt-3<i and C-2p orbitals. 



B. Adsorption of H2 to a free Pt atom 

We now consider the interaction between H2 molecules 
and a free Pt atom. Optimized binding structures are 
shown in Fig. Upon approaching to a free Pt atom, 
single H2 molecule dissociates and form PtH2 with Pt-H 
bond distance dpt_#=1.51 A and cLh-h = 2.08 A. The 
total energy Ey relative to the energies of free Pt and H 



(a) (8,0) SWNT (b) 



d H =0.95A 




E=-S90.<>46 cV 



E=-590.940 eV 



(c) GRAPHITE (d) 



,<!,= 1.87A 




E=-455.22 eV 



FIG. 5: Optimized adsorption geometry of H2 molecule on a 
single Pt atom, (a) Pt atom is adsorbed near the H-site of 
(8,0) SWNT (side and top view) (b) Pt at the A-site (bridge 
position) of (8,0) SWNT (side and top view) (c) Pt atom is 
adsorbed near the H-site of the graphite surface (d) Pt atom 
at the A-site of graphite. Et is the total energy relative to 
the constituent C, Pt and H atoms. 



atoms is calculated -9.88 eV. The binding energy relative 
to H2 molecule and free Pt atom is Eb=3.09 eV. As for 
the adsorption of H2 to PtH2 , there are several minima 
on the Born-Oppenheimer surface: The first minimum 
corresponds to a configuration in Fig. 0Jb) where PtH2 
preserves the dissociated configuration while second H2 
is molecularly adsorbed. Even if H-H interaction of ad- 
sorbed H2 is weakened and hence the H-H distance has in- 
creased to 0.87A, we identify it as molecular adsorption. 
We denote this configuration as PtH2+H2. The adsorbed 
H2 molecule is perpendicular to the plane of PtH2 . The 
binding energy of the second H2 to PtH2 is calculated to 
be Ef,=0.87 eV. Under these circumstances the average 
binding energy of each H2 is 1.98 eV. In the second config- 
uration, identified as Pt-2H2 as shown in Fig.^Jc), both 
H2 are molecularly adsorbed. As H-H molecular bonds 
are weakened, cLh-h is increased to 0.96A and all Pt-H 
bonds have uniform length with dp t -H =1.65A. Here ad- 
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sorbed H2 molecules are perpendicular. The binding en- 
ergy of each molecules is calculated to be 1.94 eV slightly 
less then the average binding energy in PtH+H2 configu- 
ration. The configuration PtH2+2H2 shown in Fig.^Jd) 
involves the adsorption of three H2 molecules; one disso- 
ciatively, remaining two are molccularly adsorbed. Here 
Pt-ti2 planes of two molecularly adsorbed H2 are per- 
pendicular. The binding energy of the third H2 rela- 
tive to the energy of PtH2+H2 in Fig. EJb) is found to 
be 0.81 eV. Accordingly, the average binding energy of 
each H2 is 1.6 eV relative to free H2 and free Pt atom. 
Fig. 0Je) compares four distinct configurations related 
with the adsorption of two molecules on the same free Pt 
atom. It appears that these configurations correspond to 
local minima on the Born-Oppenheimer surface and the 
configuration in Fig.^Jb) appears to have lowest energy. 



C. Adsorption of H2 to a Pt atom on SWNT 

We deduced two configurations for the adsorption 
of a single molecule to a Pt atom adsorbed on the 
(8,0) SWNT as described in Fig. While these two 
chemisorption configuration look dramatically different, 
their total energies differ only by 6 meV (that is smaller 
than the accuracy range of DFT). In the configuration 
described in Fig. OJa) H2 is dissociatively adsorbed with 
binding energy E&=1.1 eV relative to free H2 and Pt 
adsorbed on SWNT, i.e. E T [H 2 ] and E T [Pt+SWNTl, 
respectively. The H-H and Pt-H distances are 1.86 A 
and 2.3 A, respectively. Whereas in the configuration in 
Fig. OJ b) H2 is molecularly adsorbed with a significantly 
weakened H-H bond. H2 approaching from different di- 
rections and angles results in a chemisorption state with 
binding energy E b =1.1 eV and Pt-H distance 1.7 A. The 
length of H-H bond has increased from 0.75 A to 0.95 A 
upon adsorption 27 . Whereas, this configuration appears 
to be less energetic than the former one, the energy dif- 
ference between them is only 6 meV. Notably, while in 
the first configuration leading to dissociative adsorption 
Pt atom is located near hollow H-site, in the "molecular" 
adsorption of H2 Pt atom is adsorbed at the A-site. As 
compared to the average binding energy calculated for 
H2 adsorption to a free Pt atom, the molecular adsorp- 
tion of H2 to Pt atom adsorbed on SWNT is relatively 
weaker due to the Pt-SWNT bond. 

Adsorption of H2 to a single Pt atom attached to the 
surface of graphite is of interest in order to reveal how 
the binding energy and binding configuration of H2 de- 
pends on the radius of SWNT. We considered two config- 
uration, namely single Pt atom is adsorbed near hollow 
H-site, as shown in Fig. |3Jc) and Pt at the A-site as 
shown in Fig. Efd). For both location of Pt atom on the 
graphite surface, H2 molecule approaching the adsorbed 
Pt atom is dissociated and eventually Pt atom formed 
two Pt-H bonds with individual H atoms. In this case the 
binding with the graphite surface is weaker than that on 
SWNT, and thus dpt-c is increased to 2.4A. Relatively 
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FIG. 6: Optimized structure of H2 molecules adsorbed to 
the Pt atom on the SWNT. (a) One H 2 adsorbed to PtH 2 . 
The inset show the regions of charge depletion (Ap < 0) and 
charge accumulation (Ap > 0) as a result of the bonding be- 
tween SWNT and PtH2+H2 in (a), (b) Another local minima 
where two H2 is molecularly adsorbed to the Pt atom, (c) One 
H2 is chemisorbed two H2 are weakly bound, (d) four H2. 



weaker interaction between Pt and graphite surface al- 
lows stronger interaction between H2 and Pt, as in the 
case of free Pt atom, and hence leads to the dissociation 
of the molecule. In view of two limiting case in Fig. [5] one 
can expect that dissociation of H2 may occur for SWNTs 
having larger R. 




(d) (e) (f) 




FIG. 7: Dissociative adsorption of single H2 on a small Pt cluster adsorbed on SWNT (a). One H2 is approaching two adjacent 
Pt atoms adsorbed on SWNT. (b) Optimized geometry after dissociative adsorption of H2. (c) Variation of total energy with 
distance z. Small circles corresponds to Et for unrelaxed H2 and unrelaxed SWNT. Dotted curve corresponds to Et of the 
geometry relaxed at certain distances, (d), (e) and (f) same as (a), (b) and (c) except the Pt cluster consists of 3 Pt atoms. 2 
is the distance from the surface of SWNT. Variation of Et(z) is amplified by inset. 



To answer the questions of how many H2 molecule can 
be attached to a single Pt atom we perform a systematic 
study outlined in Fig. First we let a second H2 ap- 
proaches to Pt atom on SWNT that already has one H2 
molecule attached as in Fig. GJb) . The final optimized 
geometry of Pt and adsorbed H2 molecules in Fig. EJa) 
is similar to the configuration PtH2+H2 described in 
Fig- Elk)- This situation is explained by the weakening 
Pt-C bond between PtH 2 and SWNT shown in Fig.^a) 
which is followed by the increase of dp t ~c- First H2, 
which was initially chemisorbed to Pt as a molecule has 
dissociated upon the molecular adsorption of the second 
H2 . The dissociation of H2 is an indirect process and is 
mediated by the weakening of the Pt-C bonding through 
the molecular adsorption of second H2. The difference 
charge density Ap(r) = p T (r) - pswntO) - pptH 2 +H 2 (r) 
calculated from the difference of total charge density 
pr{v) of SWNT+PtH 2 +H 2 in Fig. EJa) and those of 
SWNT and PtH 2 +H2 indicates that while the charge of 
Pt-d xy and C-p XtV orbitals are depopulated, the Ytd z 2 
and C-p z orbitals become populated to form Pt-SWNT 
bond. This result suggests that in a reverse situation the 
weakening of the Pt-SWNT bond would lead to the trans- 
fer of charge from Pt-C bond to Pt-H bonds resulting in 
increased population of d^-orbitals in favor of the dis- 
sociation. Interestingly, exactly the same configuration 
has been obtained even when two H2 molecule approach 
concomitantly the bare Pt adsorbed on SWNT. Another 



configuration related with two molecularly adsorbed H2 
is shown in Fig. EJb), which appears to be a local min- 
ima on the Born-Oppenheimer surface and ~0.4 eV less 
energetics than that in Fig. Efa) . The configuration of 
Pt+2H2 has remained even after SWNT is removed. 

In Fig. Etc), two H2 approaching from both sites of 
PtH 2 already adsorbed on SWNT have been attached by 
weak physisorption bonds resulting PtH2+2H 2 . Their 
distances to Pt atom are relatively larger (2.lA and 3.2A) 
than that occurred for molecularly chemisorption of H2. 
The latter Pt-H 2 distance is too long and the binding en- 
ergy is ~ 20 meV; the binding energy can increase slightly 
by the VdW interaction but the adsorbed molecule can 
desorb and escape from Pt at high-temperature. Note 
that due to weak interaction between PtH 2 and both H 2 
molecules, the Pt-C bond becomes stronger, and conse- 
quently H-H distance of PtH2 is closed to be associated 
in H 2 . The attempts to attach more than three molecules 
to the Pt atom have failed. For example, as shown 
in Fig. EJd), from four H2 brought at the close prox- 
imity of Pt atom, only three were attached (one disso- 
ciatively chemisorbed, one molecularly chemisorbed, one 
physisorbed and the forth escaped). At the same time the 
Pt-SWNT bond has weakened and hence dpt-c distance 
has increased to 2.3A. We define the interaction energy 
between Pt+nH 2 , where n=2 (one of H 2 is dissociated), 
3 and 4 (one of H 2 is dissociated) and SWNT in Fig.[BJa- 
c) as E[Pt-nH 2 ]+E T [SWNT]-E T [Pt+nH2+SWNT]. Here 



7 



the total energies are calculated using the same atomic 
structures in Fig. [fja-c). Calculated interaction energies 
for each case are 0.68 eV, 1.88 eV and 0.78 eV, respec- 
tively. Using the similar procedure we also calculated the 
interaction energy between Pt+H2 (where H2 is molec- 
ularly adsorbed) and SWNT in Fig. Eta) to be 1.93 eV. 
Note that the variation of these energies with structure 
and Pt-C distances confirm the above arguments related 
with the dissociation of one of H2 followed by the weak- 
ening of the bond between Pt and SWNT. 



D. Adsorption of H2 to a small Pt cluster on 
SWNT 

As shown in Fig. Ufa), the situation is different in the 
case of interaction between H 2 and a small Pt cluster 
(consisting of a few Pt atoms adsorbed at close proxim- 
ity). As H2 approaches two Pt atoms on SWNT it starts 
to dissociate at a distance ~ 3. 9 A from the surface of 
SWNT. The optimized configuration is shown in Fig.[7Jb) 
where H-H molecular bond is broken and each adsorbed 
Pt atom formed a Pt-H bonds with dp t -H — l.mk. Upon 
chemisorption dp t -c increased from 2.15A to 2.36A. The 
dissociation process schematically shown in Fig. He) by 
plotting the variation of total energy Et as a function 
of z for two different cases. The curve by small circles 
corresponds to the total energy of SWNT+Pt 2 and H2 
calculated for different H2-tube distance z by keeping the 
atomic configuration at z — > 00 frozen for all z. The con- 
tinuous curve is obtained by relaxing the atomic configu- 
ration at certain values of z. We see that for z < 4.2A Et 
starts to lower upon the onset of dissociation. We note 
very small barrier at about z ~ 4.5A. Upon overcoming 
this energy barrier, the process is exothermic with an en- 
ergy gain of ~ 1.2 eV. As described in Fig. Efd-f), the 
adsorption of single H2 on a Pt cluster consisting of three 
Pt atoms also results in dissociation of the molecule. As 
the size of cluster increased by inclusion of the third Pt 
atom, the small potential barrier at z ~ 4. 5 A is further 
lowered, the binding energy increased to 1.5 eV. Also one 
of the Pt atoms which binds both H atom is detaches 
from the SWNT surface. This situation confirms that 
Pt-SWNT bonds are weakened upon the (molecular or 
dissociative) adsorption of H 2 to Pt. 

The interaction between Pd atoms adsorbed on SWNT 
and H 2 molecule is somehow similar to that with Pt atom. 
However, the latter case leads relatively less strong inter- 
action and smaller binding energies. For example, the 
interaction between H 2 and a single Pd atom adsorbed 
on SWNT results in a binding between chemisorption 
and physisorption with a binding energy of 0.6 eV. In 
this case, while the H-H bond length is stretched a little 
from the normal value 0.7 A to 0.8 A, the C-Pd bond is 
stretched from 2.1 A to 2.2 A. Small changes after the 
adsorption of H2 are manifestations of relatively weak 
H 2 -Pd interaction. In contrast to adsorbed two Pt atom 
in Fig. E^a) breaking the H2 molecule, two adsorbed Pd 



atoms give rise to chemisorption of molecule with more 
stretched H-H bonds. 



V. CONCLUSIONS 



In this work we presented a detailed analysis of the in- 
teraction between hydrogen molecule and a SWNT. We 
found that the binding energy between H 2 and outer sur- 
face of a bare SWNT is very weak and the physisorption 
bond can easily be broken. We showed that the binding 
of H2 to the outer surface cannot be enhanced by ap- 
plying radial deformation to increase curvature effects at 
the site facing H2 molecule. In contrast, the interaction 
between the inner surface of (8,0) tube and H2 is repul- 
sive which can prevent molecules from entering inside the 
tube. The repulsive interaction may turn to be attractive 
for large tube radius. To promote H 2 uptake on SWNT 
surface we considered functionalizcd tubes through ad- 
sorption of foreign atoms. The binding energy of H2 on 
SWNT surface did not increase by the coadsorption of Li. 
However, the binding energy increased if H 2 is directly 
attached to adsorbed Li; yet the nature of the bonding 
remained physisorption. 

The situation with Pt atom, which can make strong 
chemisorption bonds with the outer surface of SWNT is 
found to be interesting from the point of view of H 2 stor- 
age. We showed that H 2 molecule can form chemisorp- 
tion bonds with free Pt as well as Pt adsorbed on SWNT. 
Single H2 adsorbed on a free Pt atom dissociates and 
forms two strong Pt-H bonds. On the other hand, while 
single H 2 molecule is molecularly chemisorbed to a single 
Pt atom at the A-site of SWNT surface, it can dissoci- 
ate if Pt atom adsorbed near the hollow site. Even the 
molecular adsorption of single H2 can turn dissociative 
if a second H2 is molecularly adsorbed to the same Pt 
atom. The dissociative adsorption is mediated by the 
weakening of Pt-C bonds either due to a specific location 
of Pt on SWNT or due to the second H2 molecularly ad- 
sorbed to Pt. Dissociative adsorption of single H2 to a 
single Pt atom on the graphite surface suggests that the 
dissociation of H 2 is favored on SWNTs having large ra- 
dius. Our analysis suggests that single Pt adsorbed on 
SWNT can bind up to two H 2 molecules with significant 
binding energy in the chemisorption range. Beyond two 
adsorbed H2, additional molecules form weak physisorp- 
tion bonds with Pt. Single Pd atom adsorbed on SWNT 
exhibits similar effects but in relatively weaker manner as 
compared to that of Pt. Interesting interaction between 
H 2 and Pt and resulting bonding mechanisms justifies 
similar investigations of SWNTs functionalized by other 
transition elements (in particular Ti, Ni, Cr, V) for a 
higher wp H2 storage. 
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